Evidence for Kosterlitz-Thouless type orientational ordering of 
CFaBr monolayers physisorbed on graphite 
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Monolayers of the halomethane CFsBr adsorbed on graphite have been investigated by x-ray diflrac- 
tion. The layers crystallize in a commensurate triangular lattice. On cooling they approach a three- 
sublattice antiferroelectric pattern of the in-plane components of the dipole moments. The ordering 
is not consistent with a conventional phase transition, but points to Kosterlitz-Thouless behavior. 
It is argued that the transition is described by a 6-state clock model on a triangular lattice with 
antiferromagnetic nearest neighbor interactions which is studied with Monte-Carlo simulations. A 
finite-size scaling analysis shows that the ordering transition is indeed in the KT universality class. 



PACS numbers: 61.10.Nz, 64.60. Cn, 68.35.Rh 

Two-dimensional systems cannot show long-range or- 
der by breaking a continuous symmetry at any finite tem- 
perature T. Nevertheless Kosterlitz and Thouless (KT) 
demonstrated that the 2D XY-ferromagnet with nearest 
neighbor interaction and equivalently the planar rotator 
exhibit a phase transition from an ordered phase with an 
algebraic decay of the spin correlations to a disordered 
phase via the unbinding of vortex-antivortex pairs Q . 

The KT-concept has been applied successfully to the 
onset of superfluidity |^,^ and of superconductivity in 
thin films and two-dimensional melting An applica- 
tion to layered magnetic systems with large ratios of the 
intra- and the inter-layer exchange appears to be more 
direct, but here the KT transition is masked by the tran- 
sition to 3D behavior . In the present article we present 
x-ray diffraction data on CFsBr monolayers physisorbed 
on graphite and present evidence that the orientational 
pseudospin ordering is of the KT-type. To our knowl- 
edge the present study reports the first observation of a 
KT-transition in a (quasi)-magnetic system. 

Physisorbed monolayers are indeed good approxima- 
tions to two-dimensional (2D) systems, as evidenced by 
the fact that phase transitions occurring in such layers 
can be often described by the 2D versions of elementary 
models of statistical mechanics such as the Ising or the 
3-state Potts model [0-§. 

The halomethane CFsBr is a prolate C^y molecule with 
a dipole moment of about 0.5 D. Monolayers have been 
adsorbed on exfoliated graphite and investigated by x- 
ray powder diffraction. The coverage p, temperature T 
phase diagram is rather complex Here we con- 

centrate on a coverage which is representative of the ex- 
tended monolayer regime in which the monolayer lattice 
is commensurate with the graphite lattice. Diffraction 
patterns are shown in Fig. 1. The peaks of the patterns 
have been analyzed for position, intensity, and the coher- 
ence length, as has been described elsewhere The 
commensurate layer melts around 105 K p2|. In the 2D 



solid state, the diffraction pattern is dominated by the 
principal Bragg reflection of the commensurate, triangu- 
lar 2 x 2 lattice at Qp= 1.475 A~-^. Just below the melt- 
ing temperature, this peak is in fact the only one visible. 
Down to lowest T, the peak position is independent of T, 
as imposed by the commensurability with the substrate. 
The peak width yields a coherence length of about 180 A 
throughout the 2D solid state which compares well with 
the average lateral size of the crystallites of exfoliated 
graphite usually quoted in the literature. The 2x2 struc- 
ture has been also observed for monolayers of CF4 , 
CF3CI (l|], C2F6 Hi on graphite. These molecules and 
CFsBr have the threefold symmetry and the CF3 group 
in common. Hence it is plausible to assume that CFaBr 
stands on the substrate, presumably with the F3 tripod 
down. An isolated CFsBr would prefer to lie flat on the 
substrate, but definitely the 2x2 mesh is too tight to 
accommodate the molecules in this orientation. On the 
other hand steric considerations suggest that the lateral 
2x2 packing of perpendicular molecules can tolerate tilt 
angles of the molecular axis up to 30° with respect to the 
substrate normal. A tilt is equivalent to a non-zero in- 
plane component of the dipole moment. We regard this 
component as pseudospin Si . In this sense the 2x2 state 
is disordered with a zero time average of every Si, and is 
stabilized at higher T by a libration and/or a precession 
of the molecular axis about the substrate normal. 

On cooling, two additional features develop in the 
diffraction pattern which finally, below 40 K, can be re- 
garded as Bragg peaks at Q=0.851 and 2.253 A~^. 
They are easily identified as superlattice reflections 
brought forth by a tripling of the 2D mesh leading to 
a 2V3 X 2-^/3 commensurate hexagonal lattice. The (hk) 
indices of Fig. 1 refer to this superlattice. The (20) reflec- 
tion, which should occur at 1.703 A~^, is absent. Below 
40 K, the coherence length derived from the width of the 
superlattice reflections is again 180 A and is thus imposed 
by the substrate. The global structure of the low temper- 
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FIG. 1. Diffraction patterns of solid monolayers of CFsBr 
on graphite (p = 1.08). The (hk) indices refer to the 
2\/3 X 2\/3 superlattice. The blind region is due to the strong 
(002) reflection of the graphite substrate. 



ature state must be consistent with the hexagonal plane 
group p3 with three molecules in the supercell. Regard- 
ing the molecules as rigid, there are five positional param- 
eters, namely the translational coordinates x,y and the 
three Eulerian angles of one molecule. The low number 
of reflections does not allow a rigorous determination of 
these parameters, nevertheless one can make some semi- 
quantitative statements. A perpendicular orientation of 
the molecules in combination with a staggered azimuthal 
angle of the CF3 group about the substrate normal is 
not sufficient to explain the superlattice intensities. The 
cell tripling must involve tilts of the molecular axes away 
from the perpendicular direction, presumably in combi- 
nation with slight displacements of the molecule centers 
out of the original positions. A satisfactory agreement 
of measured and calculated intensities (including the ac- 
cidental absence of the (20) reflection) is obtained with 
tilts of about 20° pointing toward a next neighbor such 
that the pseudospins form an antiferro- "magnetic" 120° 
Neel-pattern on a 2D triangular lattice. 

In case of a conventional order-disorder transition, the 
integrated intensities of the superlattice reflections would 
be a measure of the order parameter and should hence 
vanish when approaching the critical temperature from 
below. (See ref. |^ for such a transition for C2F5CI on 
graphite). This is however not what is observed in the 
experiment on CFsBr monolayers. On heating the super- 
lattice peaks broaden while keeping approximately their 
integral intensity until they are finally lost in the back- 
ground. The temperature evolution of the (21) reflection 
is shown in Fig. 2 together with fits of a theoretical pro- 
file. The T-dependence of the coherence length ^ as de- 
termined from the intrinsic width of this peak is shown 
in Fig. 3. The solid line of this figure is a fit of the KT- 
expression 
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FIG. 2. The temperature evolution of the (21) superlattice 
reflection. 



to the data for T > 40 K. The fit parameters are A = 
9±2 A, B = 1.5±0.4, Tk = 30±3 K. Note that the value 
of A is reasonably close to the lattice parameter of the 
2D mesh. Thus ^ is expected to diverge at a KT-critical 
temperature Tk of about 30 K, but the growth of the 
correlated regions is interrupted when ^ reaches the size 
of the graphite crystallites. This happens at about 40 K. 

Given the scatter of the data points of Fig. 3, it is clear 
that fits with power laws ^ oc {T/Tc — 1)^'' are likewise 
possible. We do think, however, that an interpretation of 
the intensity between 40 K and 55 K in terms of conven- 
tional critical diffuse scattering, T > Tc, is not meaning- 
ful, because of the absence of any critical T-dependence 
of the integrated intensities, and because the diffuse scat- 
tering would be presumably far below the detection limit 
in x-ray or neutron diffraction studies on adsorbed mono- 
layers on exfoliated graphite. At least we do not know of 
any successful experiment of this type. 

The heat capacity measurement also shows an indica- 
tion of the phase transition |12|. Apart from the melt- 



ing anomaly, the heat capacity of commensurate CFaBr 



2 



^ 200 I I I i I I ■ I I I 




CD - ' - 

"5 Q L_! 1 I I I 1 I I I I I I 1 I 

o 30 40 50 60 

temperature T (K) 

FIG. 3. The temperature dependence of the coherence 
length of the antiferroelectric correlations, as derived from 
the intrinsic width of the (21) reflection. 

monolayers shows an anomaly extending from 35 K to 
40 K which is just the T-range where ^ saturates. How- 
ever it does not allow a quantitative analysis of the crit- 
ical scaling behavior since the anomaly is very weak. 

Do CFsBr monolayers meet the theoretical require- 
ments for a KT-transition? Clearly the pseudospin cor- 
relations are bound to a plane, thus the system is 2D 
with respect to the relevant degrees of freedom at the 
phase transition. The pseudospin is presumably not a 
strictly isotropic planar rotator but experiences a crystal 
field which breaks the continuous azimuthal symmetry. 
It has been shown however that — at least for the planar 
XY ferromagnet — the KT-behavior is stable with re- 
spect to small crystal fields of 6-fold symmetry jl^ . The 
site symmetry of the monolayers is indeed 6-fold, but the 
ordered structure approached is not ferro- but antiferro- 
electric. Antiferro-type ordering on a triangular lattice 
is affected by frustration with the consequence that the 
helicity, a discrete 2-fold symmetry, enters into the prob- 
lem. The order-disorder transition is then described by a 
confluence of the Ising and KT universality classes [ |l7|Jl8| 
with an Ising type anomaly of the specific heat, but the 
spin correlations of the disordered phase still follow the 
KT-form of ^ (T) as we will demonstrate now. 

From the above considerations follows that the pseudo- 
spins formed by the CFsBr dipoles are arranged on 
a triangular lattice, interact antiferromagnetically with 
nearest neighbors and have 6 preferred orientations. 
Hence the model that captures the universal features of 
the orientational ordering transition should be the anti- 
ferromagnetic six-state clock model on a 2D N = L^x Ly 
triangular lattice (see Fig. The six-state clock spin 
is a planar spin pointing toward discrete six directions; 
S = (cos6',sin6') with 61 = ^ (n = 0, 1, . . . , 5). The 
Hamiltonian of this model reads 

n = 2j^cos{e^-ej) , (2) 




FIG. 4. Left: The spatial arrangement of the CFsBr 
molecules on the underlying hexagonal graphite surface in 
the 2x2 structure that is present below 105K. The resulting 
trigonal lattice is also shown. Right: A triangular lattice of 
size 12 X 3 with the 6-state clock spins indicated as arrows 
representing the 6 preferred directions of the CF3B dipolar 
moments. The arrow configuration shown is one of the 12 
ground states by which the lattice sites are separated into 
three A, B, and C sub-lattices. The helicity for each ele- 
mentary triangle corresponding to this ground state is also 
indicated by + and — . 

where the sum is over nearest neighbor site pairs 
and J > is the antiferromagnetic coupling strength. 
The overall factor 2 is introduced for a computational 
convenience. It is a limiting case of an antiferromagnetic 
XY model with infinite anisotropy field cos p6 (p — 6). 

The antiferromagnetic coupling on a triangular lattice 
induces a frustration. As a result the system has 12-fold 
degenerate ground states with three-sublattice structure. 
Global spin rotations by multiples of 27r/6 and also Ising 
like transformation that inverts all helicities formed by 
the spins on elementary triangles (see Fig. 4) transforms 
one ground state into another. Therefore the model pos- 
sesses a Cg (6-fold clock) x Z2 (Ising) symmetry and we 
expect that there are two phase transitions associated 
with the spontaneous breaking of the Cg and Z2 symme- 
try. 

We performed extensive Monte-Carlo simulations of 
the model (|) using conventional single spin-flip dynam- 
ics with up to 10'* Monte-Carlo steps per spin and lattice 
sizes up to 384 X 192 and studied its phase transition with 
a thorough finite size scaling analysis. Here we focus on 
the breaking of the Cg symmetry that is indicated by the 
x-ray diffraction experiments reported above, the details 
of the breaking of the Z2 symmetry (which happens at 
a slightly higher temperature and turns out to be in the 
Ising universality class) will be reported elsewhere po| . 

The order parameter for the Cg symmetry is the sub- 
lattice magnetization 

rriA = ^^exp{i9,) , (3) 

where the sum runs over only sites in the A sublattice 
{niB and mc can be defined similarly). Analogously the 
spatial spin correlation function is defined as 
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FIG. 5. Left: The Binder parameter Bm (P) versus tem- 
perature for different system sizes. The intersection point of 
the different curves locates the critical temperature indicated 
by the broken line at Tc^ = 1.035±0.0005. Middle: The cor- 
relation length as a function of T — Tcg (with Tcg = 1.035) 
in a log-log plot. Upwards curvature of the numerical data 
indicate already a faster-than-algebraic divergence. The full 
line is a least square-fit to the KT-form (0). Right: The 
magnetic susceptibility with a fit to the expected KT-form. 
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and the susceptibility is given by Xm = N{m\). 

The transition temperature is estimated using the di- 
mensionless ratio of moments (or Binder parameter) 
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which is shown in Fig. 5 and yields the critical temper- 
ature — 1.035 ± 0.0005. We extract the correlation 
length from the correlation function C{r) defined in (^) 
via 



(6) 



Our data are shown in Fig. 5 in a log-log plot of £m(T) 
versus T — Tcg with Tcg from above. They fit nicely 
to the KT-form (0) with A = 2.54 and B = 0.193. 
Note that the non-universal number, B, differs signif- 
icantly from the one for the experimental data, which 
is not unusual for microscopically different systems in 
the KT universality class (see e.g. |2^). The susceptibil- 
ity, also shown in Fig. 5, follows the expected KT-form 
Xm = A'exp{B'{T - TcJ"^/^), too. We also checked 
other quantities (like the decay of Cm(r) at Tcg) and 
found everything to be consistent with a KT-scenario. 

In summary, our measurements show that CFsBr 
molecules physisorbed on graphite are arranged in a com- 
mensurate triangular lattice. On cooling, a Neel-type 
120° pattern is approached. The growth of the criti- 
cal correlations has been followed over a wide tempera- 
ture range. It finally stops when the correlation length 
reaches the size of the substrate crystallites. The behav- 
ior is consistent with the KT-scenario. An antiferromag- 
netic 6-state clock model on a triangular lattice, which 



we studied numerically, describes the universal features 
of this ordering transition well. 
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